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ABSTRACT 

We present a simple model for the relationship between quasars, galaxies, and dark matter halos from < 
z < 5. In the model, black hole (BH) mass is linearly related to galaxy mass, and galaxies are connected to dark 
matter halos via empirically constrained relations. A simple "scattered" light bulb model for quasars is adopted, 
wherein BHs shine at a fixed fraction of the Eddington luminosity during accretion episodes, and Eddington 
ratios are drawn from a lognormal distribution that is redshift-independent. This model has two free, physically 
meaningful parameters at each redshift, the normalization of the MaH-A^gai relation and the quasar duty cycle; 
these parameters are fit to the observed quasar luminosity function (LF) over the interval 0.5 < z < 4.75. This 
simple model provides an excellent fit to the LF at all epochs, and also successfully predicts the observed 
projected two-point correlation of quasars from 0.5 < z < 2.5. It is significant that a single quasar duty cycle 
at each redshift is capable of reproducing the extant observations. The data are therefore consistent with a 
scenario wherein quasars are equally likely to exist in galaxies, and therefore dark matter halos, over a wide 
range in masses. The knee in the quasar LF is a reflection of the knee in the stellar mass-halo mass relation. 
The normalization of the Mbr-M^^i relation increases as (1+z)^, consistent with observations, although in 
our model the normalization is degenerate with the mean Eddington ratio. Future constraints on the quasar LF 
and quasar clustering at high redshift will provide strong constraints on the model. A further implication of 
our model is that the autocorrelation function of quasars becomes a strong function of luminosity only at the 
very highest luminosities, and will be difficult to observe because such quasars are so rare. Cross-correlation 
techniques may provide useful constraints on the bias of such rare objects. The simplicity of our model allows 
for rapid generation of quasar mock catalogs from N-body simulations that match the observed luminosity 
function and clustering to high redshift. 

Subject headings: quasars: general — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Quasars are among the most luminous astrophysical ob- 
jects, and are believed to be powered by ac cretio n onto 
supermassive black holes (e.g. ISalpeterlll964l: iLvnden-Belll 
Il969i) . They have become a key e lement in our cur- 
rent paradigm of galaxy evolution (e.g., iSpringel et al.l 120051 : 
ICroton et alJl2006t iHopkins e t al.ll2008h . and essentially all 
spheroidal systems at present harbor massive black holes 
dKormendv & Richst one 1995), the masses of which are cor- 
related with many properties of their host systems. Despite 
their importance, and intense theoretical activity, a full theory 
of the coevolution of galaxies and quasar eludes us. 

The current paradigm assumes that every galaxy ini- 
tially forms in a gas-rich, rotationally-supported sys- 
tem. Once the dark matter halo grows to a critical 
scale s ome event, mos t likely a major merger (ICarlbergj 
1990^, 'Haiman & Loebl 11991 ICattaneo, Haehnelt & Re^ 
1999; Kauffmann & Haehneh 2000; Springel et al. 2001 
Hopki ns et al. 2006, 2008) or instability in a cold-stream fed 
disk (Ciotti & Ostriker 1997, 2001; Di Matteo et al. 201^, 
triggers a period of rapid, obscured star formation, the gener- 
ation of a stellar bulge and a growing black hole (BH). Even- 
tually the accreting BH becomes visible as a quasar, and soon 
after the star formation is quenched on a short timescale, per- 
haps via radiative or mechanical feedback from the BH (e.g. 
Silk & Rees' 1998; King 2003; Wvithe & Loeb 2003; Shankar 
2009; Natarajan 2012;' Alexander & Hickox 2012). Under- 
standing the details of this picture remains an active area of 



research. 

Phenomenological models for quasar demographics often 
adopt power-law relation s between quasars, galaxies, and 
dark matter halos (e.g., Efstathiou & ReesI 11988'; Carlberd 
1990; Wvithe & Loeb 2002, 2003; Haiman, Ciotti & Ostrik^ 
2004; Marulli et al. 2006; Lidz et al. 2006; iCrotonI I2009r 
Shen 2009; Booth & Schaye .2010 ). In these models, the 



duty cycle of quasars is tuned to match the observations, 
and a generic conclusion is that the duty cycle is a strong 
function of halo mass or quasar luminosity, peaking at a 
halo mass of 10'^~'^Mq. However, these previous models 
do not incorporate constraints provided by the galaxy stellar 
mass function over the interval < z < 5. And yet, a 
variety of lines of evidence suggest that the relation between 
halos and galaxies is highly non-linear, with a characteris- 
tic peak in galaxy formation efficiency at a halo mass of 
~ lO'^MfT, dvan den Bosch et al.| llOOl |We"& Ostrike^ 
12004; Mandelbaum et al. 2006; Conrov & Wechslep 
gp09; Moster et al. .2010^ .Truiillo-Gomez et aL ,201 iF 
Eehroozi, Wechsler & Conrovll2012h . The aim of this paper 
is to incorporate empirically constrained relations between 
galaxies and halos into a simple model for quasar demo- 
graphics. We will demonstrate that a model constructed to 
match the observed galaxy stellar mass function implies a 
quasar duty cycle that is independent of galaxy and halo mass 
at each redshift. This has important implications for physical 
models aimed at understanding the triggering of quasars and 
their connection to the evolution of galaxies. 
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The outline of the paper is as follows. In ^we describe 
the model, in ^ the model is compared to data, and a dis- 
cussion is presented in ^ We conclude in ^ Where neces- 
sary we adopt a ACDM cosmological model with ri„, = 0.28, 
ft^ = 0.72 and erg = 0.8. Unless the h dependence is explic- 
itly specified or parametrized, we a ssume /; = 0.7. Dark mat - 
ter halo masses are quoted as M^i dBrvan & Normanlll998l) . 
Luminosities are quoted in Watts an d magnitudes in the AB 
system, and stellar masses assume a IChabrieri (l2003l) stellar 
initial mass function. 

2. THE MODEL 

Our goal is to construct a simple model that relates galaxies, 
quasars, and dark matter halos over the redshift interval < 
z < 5. A small number of free parameters will characterize 
the model, and these parameters will be constrained against 
observations. 

We begin by specifying a dark matter halo mass function 
an d its evolution to z = 5. We adopt the fitting functions 
of lTinker etall (120081 12010') for the halo mass function and 
large-scale bias, which represent the latest fits to these pa- 
rameters from cosmological A^-body simulations. Note that 
here and throughout we consider only parent halos; satellite 
halos, also known as subhalos, or not included in the present 
study. This is a reasonable approximation at high redshift, as 
quasars inhabit highly biased halos on the steeply falling tail 
of the mass function and any satellite galaxies of the same 
mass would live in even more massive halos which are ex- 
ponentially rare. This assumption will break down at lower 
luminosities, where the satellite fraction can be expected to 
rise. This assumption will also fail to account for the small- 
scale clustering of quasars, in particular the clustering within 
the halo scale of < 1 Mpc. When we compare to clustering 
measurements in i l3.2l we will therefore restrict our compari- 
son to R > 1 Mpc, which is where most of the data lie. 

Next, we adopt empirically constrained relations be- 
tween galaxy stellar mass and da rk matter halo mass over 
the in terval < z < 5 from Beh roozi. Wechsler & Conrovl 
(12012b . These relations were constrained against the ob- 
served galaxy stellar mass function and galaxy star forma- 
tion rates from z = to z = 8. The resulting stellar mass- 
halo mass (Mgai-M/,) relation agrees with results obtained 
from other techniques, including abundance matching, halo 
occupation models, satellite kinematics, and gravitational 
lensing (see Behroozi, Wechsler & Conroy 2012). We also 
adopt an amount of scatter between galaxy mass and halo 
mass as a function of redsh i ft imp lied by the model of 
iBehroozi. Wechsler & Conrovl (120121) . This scatter increases 
from « 0.2 dex at z = 0.5 to « 0.4 dex at z = 4. 

Galaxies are assigned BHs via the following equation: 

where Mgai and Mbh are the stellar mass of the galaxy and 
mass of the BH, respectively. The available data at z ^ 
is consistent with a linear relation between Mgai and Mbh, 
(i.e. /3 = I) which is what we adopt herein, with a normal- 
ization constant of a « -3.1 (Haring & Rix 2004). The scal- 
ing with redshift is motivated by obser vations (IMcLure et alj 
l2006tlTargett. Dunlop & McLurel2012h . but since we fit for a 
at each redshift, any deviation from (1 +z)^ will be absorbed 
in the redshift-dependence of the parameter a. We adopt a 
scatter in this relation of 0.3 dex, independent of mass, con- 



sistent with the observ ed scatter in the local Mbh - c relation 
dTremaine et al.ll2002h . 

We have chosen to relate Mbh to the total stellar mass of 
the galaxy, rather than specifically to the bulge component. 
Obviously for bulge-dominated galaxies the distinction is ir- 
relevant, but the differences can grow as we include galaxies 
with a large disk component. Assuming that bulge properties 
are the dominant factor in determining Mbh, a more refined 
model would include the evolution and mass-dependence of 
the bulge-to-total ratio. However for now we neglect this dis- 
tinction. We do find that our results are relatively robust to 
modest changes in the slope of the Mbh - Mgai relation (see 
© — and any overall normalization change can be absorbed 
into our parameter a — so there are reasons to believe a more 
complex model would achieve a similar level of success in 
fitting the observations. 

In addition to the strong observed correlation between Mg^i 
and Mbh, there are well-known correlations between Mbh and 
other parameters of the galaxy includ ing the velocity disper- 
sion, a, and galaxy size, R^. In fact. iHopkins et all (l2007bl) 
argued for the existence of a BH fundamental plane (relating 
Mbh, f, and R^) that has smaller scatter than any other rela- 
tionship between Mbh and a single galaxy property. Another 
option would therefore have b een to connect BHs to galaxies 
via a, as for example done by ICrotonI (120091) . or via the BH 
fundamental plane. We choose to use Mgai herein because this 
quantity is readily available for galaxies to z = 5, and because 
the redshift-dependent connection between galaxies and ha- 
los is presently available for galaxy stellar masses, but not for 
galaxy velocity dispersions. 

The BH mass is converted to a bolometric quasar luminos- 
ity using an Eddington ratio, L/L^n = 77, that is independent 
of redshift: 

Lq = 33x\Q\^Lq. (2) 

Mq 

We draw 77 from a lognormal distribution with mean of 77 = 0. 1 
and a dispersion of 0.3 dex, in agreeme nt with observations 
dKollmeieretal .1120061; IShen et al.ll2008l) . In our model the 
value of the Eddington ratio is degenerate with the normaliza- 
tion of the Mbh - Mgai relation and any intrinsic width in the 
Eddington ratio distribution is degenerate with scatter in the 
Mbh -Mgai relation. 

In order to compare to observations, we must translate 
Lq into magnitudes in a given filter We adopt the rela- 
tion between bolometric luminosity and /- band mag nitude (k- 
corrected to z = 2) using the relation from IShen et al. (2009): 

Mi(z = 2)= 72.5-2.5 logLg (3) 

= -5.26-2.5 log(?7MBH) (4) 
= -30.3-2.5 (log T] + a)-5 log( 1 + z) 

-2.5/3 log (Mgai/lO'^Mo) , (5) 

where Lq is in Watts and Mbh is in solar masses. The last two 
relations follow directly from Equations [T] and |2l we include 
them here to make explicit the connection between Mgai and 
observed quasar magnitude, and also to emphasize the fact 
that rj and a are perfectly degenerate in our model. There is 
scatter in Lq at fixed Mgai which arises from a combination of 
scatter in Mbh -Mgai and Lq -Mbh- Assuming 0.3 dex scatter 
in each relation and adding them in quadrature we obtain a 
0.42dex scatter in Lq at fixed Mgai. 

There are two free parameters in this model at each redshift: 
the normaUzation of the Mbh -Mgai relation, specified by a. 



Quasar Demographics 



3 



logM„,, (M,J 




logM,(M,J 



Figure 1. Summary of the model relations at z = 2. The quasar LF deter- 
mines the abundance (see the points on the curve, which label space densities 
in units of log Mpc"') of quasars at a given luminosity (right vertical axis) 
or BH mass (left vertical axis). For an assumed lifetime, tQ, this maps to an 
abundance of galaxies and the stellar mass function provides the appropri- 
ate galaxy stellar mass (upper horizontal axis). The empirically constrained 
M„iii — Ml, relations from Behroozi, Wechsler & Conroy 12012) allow us to 
map this into a halo mass (lower horizontal axis). The curve shown is at 
z = 2, though the general behavior is similar at other redshifts with a steep 
low-mass slope and a shallower high mass slope (see Figure [8). Note the 
lower horizontal axis determines the clustering amplitude at fixed redshift 
while the left vertical axis determines the quasar luminosity. 

and the quasar duty cycle, /on- An important, and novel fea- 
ture of this model is that we adopt a constant duty cycle, inde- 
pendent of luminosity, Mbh or M/,. Sometimes the duty cycle 
is recast into a "lifetime" using the Hubble time: tg = fontn- 
As we will demonstrate in the following section, both of these 
parameters are highly constrained by the observed quasar LF. 

The resulting relations between galaxies, halos, and quasars 
are illustrated in Figure [T] These relations represent the best- 
fit model constrained by the quasar LF at z = 2 (see i l3.lt . The 
quasar LF allows us to relate luminosity to number density. 
For an assumed duty cycle we then have the abundance of 
BHs of that mass. Similarly the stellar mass function maps 
galaxy mass to abundance. Thus at fixed duty cycle we obtain 
a tight constraint on MBH-A^eai- As the stellar mass function 
and quasar LF contain significant curvature only one combi- 
nation of normalization and duty-cycle provides a good fit to 
the data for a range of luminosities (unless we allow signifi- 
cant variation in the lifetime as a function of luminosity). 

Figure |2] shows how the predicted quasar luminosity func- 
tion in our model (at z = 2) depends upon several parame- 
ters in the model. The amount of scatter in the Lg-Mgai 
relation is important for the shape at high luminosity, and 
indeed the abundance of luminous quasars provides a lower 
limit on the scatter for any model which places quasars in 
halos on the exponentially falling part of the mass func- 
tion. We see that a model with no scatter in the Mbh ~ 
Mgai relation predicts drastically fewer bright quasars than 
a model including scatter (see a lso White, Martini & Cohn 
2008"; 'Sha nkar. Weinberg & Shenll2010- DeGrafet aL,201k 
Trainor & Steidelll2012l for related discussion). Variations in 
the BH mass at fixed galaxy mass (a) change both the normal- 
ization and shape of the luminosity function while variation in 




M,(z=2) 



Figure 2. Variation in the predicted luminosity function of quasars at z = 2 
as a function of the parameters in our model. The dashed (red) line shows 
how the inclusion of scatter in the Mbh ~ M^.^i relation is important at the 
high mass end. with models including more scatter predicting more luminous 
quasars. Variations due to changes in the normalization of the Mbh ~Mg^i 
relation (-3.4 < a < -2.8; Equation[T) are shown by the dotted (blue) Hues, 
and we see this parameter changes both the normalization and shape of the LF 
since the galaxy stellar mass function has a particular shape. Finally the dot- 
dashed (green) line shows variation in the logarithmic slope of the Mbh "Mj^i 
relation (0.5 < /? < 1.5; Equation[T). 

the slope of the relation (/3) has a large effect on the shape of 
the LF both at low and high luminosity. 

3. COMPARISON WITH OBSERVATIONAL DATA 

3.1. The Quasar Luminosity Function 

Figure |3] shows the predictions of our model comp a red to 
a compilation of observational data from 'Wolf e t ap (I2003L 
COMBO-17; open s quare s )^|Richards et al. (20Qi SDSS; 
soHd circles). Groom etal. (2009, 2SLAQH-SDSS; open di- 
amonds), Glikm an et al. (2010, NDWFSh-DLS; stars), and 
[Masters et all (12012, GOSMOS; crosses). We have adopted 
the following transfor mation between fih ers (IWolf et al.l2003l: 
[Richards et al . 2006; Groom et al.ll2009h : 

M,(z = 2)=Mj(z = 2)-0.25 (6) 
=Mi45()-0.29 (7) 
=Mi,-0.71 (8) 

in order to convert all of the measurements to the M,(z = 2) 
system for comparison. 

The lifetime, tg, and normalization of the Mbh-M/, rela- 
tion (a in Equation [1]) have been fit to the data at each red- 
shift. The grey shaded regions mark the one sigma range of 
allowed models. In most panels the formal errors are so small 
that the grey band is buried behind the best-fit relation. The 
constraints on the parameters are so strong because the data at 
z < 4 samples luminosities both above and below the knee in 
the LF and because the formal errors on the LF are small. 

For comparison we also show the luminosity function that 
results from assuming a power-law relation between quasar 
luminosity and halo mass, as has been assumed in many 
early works (e. g. Efstathiou & Rees 1988; Garlberg 199C; 
'Wvithe & Loeb 2002, 2003; Haiman, Giotti & Ostriker 2004; 
.MarulU et al._2006;.Lidz et al..2006;.Grotoa2009;.Shea2009l; 
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Figure 3. The luminosity function of (active) quasars predicted by our model at different redshifts, as compa red to the observa tions and a simple model in 
which quasar luminosity is tied to halo, not galaxy, mass (denoted PLM for power-law model). The data are from'Wolf et al.' f2003', COMBO-17; open squares];^ 

f .chards et al. (2006, SDSS; solid circles), Croom et al. (2009, 2SLAQ+SDSS; open diamonds), Glikman et al. (2010. NDWFS+DLS; stars), and Masters et al] 
012, COSMOS; crosses). The lifetime, tQ, and the MBH~^sal normalization, a, are fit in each panel and the grey region illustrates the Icr uncertainty in 
the model prediction. Only black symbols are included in the fits; the grey symbols generally represent data of lower quality and are included for comparison 
purposes only. 

In Figure |4] we show the two derived parameters in our 
model, the quasar lifetime, fg (or, equivalently, the duty cy- 
cle), and the normalization of the Msu-Mg-^i relation, a. In 
the top panel of Figure |4] we include lines of constant duty 
cycles of 10"^ and 10"^. For reference, the Salpeter time is 
the e-folding time for a BH growing at a fraction rj of the Ed- 
dington luminosity with a radiative efficiency of e and is de- 
fined as fsaip = 4 X 1 0^(6/77) yr. It is striking how little fg varies 
from 0.5 < z < 3. The evidence for a decrease in fg at z > 3 
should be regarded as tentative, as the data used to constrain 
these parameters becomes rather uncertain and is compiled 
from heterogeneous sources. Our estimates of fg are in good 



iBooth & Schayell2010l) . This model is also characterized by 
two free parameters, the duty cycle and the normalization of 
the (power-law) relation between quasar luminosity and halo 
maso The fundamental difference between our model's pre- 
dictions and these power-law models is that we explicitly take 
into account the efficiency of galaxy formation as a function 
of mass and redshift (see Figure[TJ. The two models differ less 
significantly at higher redshifts for reasons to be discussed be- 
low. 



' The particular model we consider is Lq = -yMj^ '*, where 7 is the free 
normalization and the index, 1.4. was chosen from the power-law model of 
ICroton (2009). 
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Our Model 
McLure et al. (2006) 
z=0 




Figure 4. Upper Panel: The duty-cycle, or quasar lifetime, as a function of 
redshift. We define tQ = fontn where tn is the Hubble time at redshift z and 
/on is the probability that a BH is a luminous quasar (which is independent 
of luminosity in our model). Also shown are lines of constant /on = 10"^ and 
10"^. Lower Panel: The evolution of the normalization of the Mbh — Mg^i 
relation in our model (solid points) compared to a representative sample of 
results from the literature (see text for further discussion). 

agreement wit h quasar lifetim es inferred by other methods, as 
summarized in' Martinil (120041) . 

In the bottom panel of Figure |4] we show the evolution 
of the normalization of the Mbh - A^gai relation as inferred 
from our model. In this pa nel we also include the nor- 
malization measured at z ^ (iHaring & Rixll2004h . and es- 
timates ofjts_evolution in samples of massive galaxies to 
z ~ 2 jMcLure et al. 2006). Our model implies an increase 
in the normalization by roughly an order of magnitude by 
z ~ 2 - 3, which is in good agreement with a variety of obser- 
vational constraints to z ^ 6 (|P eng et al. 2006; Merloni et al. 
l2OT0t iDecarhet al.ll2010t llk'gett. Dunlop & McLur&.201Z) . 
Among re cent models, w e are in goo d agreement with th e 
models of iHopkins et aP (l2007al) and ICroton et all (l2006h . 
which predict roughly an order of magnitude increase in Mbh 
at Mgai ~ 1 0'° between z = a nd z = 3. In contrast, the sim- 
ul ations of ISiiacki et al.l (12007b and the semi-analytic model 
of iFanidakis et al.l ( 120121) predict almost no evolution at the 
massive end, in conflict with both our model and the data. 

Figure |5] shows the model LFs at z = 0.5 and z = 2.4. Here 
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Figure 5. Contribution to the quasar LF from quasars in different halo 
masses. The curves represent the model LF computed including halos less 
massive than the values shown in the legend (in units of logM©). The model 
LF is almost entirely insensitive to the presence or absence of quasars in halos 
more massive than 10"^Mq. 

we consider the contribution to the total LF from quasars 
in halos of different masses. Specifically, we construct 
model LFs by selecting quasars residing in halos less mas- 
sive than log(M/,/MQ) < 13.0, 13.5, and 14.0. The purpose 
of this figure is to demonstrate that massive halos contribute 
very little to the total LF. In fact, the model is almost en- 
tirely insensitive to what happens in halos more massive than 
\og{Mh/ Mq) <13.5, owing to their rarity relative to lower 
mass halos. This has important consequences for any model 
that is tuned to match the quasar LF, as we discuss in f|4] 

In Figure[3]we adopted our fiducial values for the slope and 
scatter in the MBH-Mgai and Lq-Mbh relations. We found 
that we can find equally good fits if we adjust the scatter in the 
Lg-Mgai relation between 0.3 and 0.6 dex, if we modify the 
slope of the Mbh -Mgai relation to /3 = 4/3 or 5 /3, or even if 
we change the overall normalization in the Mgai-M/, relation. 
These changes result in different best-fit values for tq and a. 
Similarly, since a and -q are perfectly degenerate in our model 
we can find very good fits if the Eddington ratio is allowed to 
vary, although large changes in 77 lead to disagreement with 
the more direct measures of the MBH-Mgai relation shown in 
Figure H] Future constraints on the Mbh - gai relation as a 
function of redshift will, in the context of our model, provide 
strong constraints on the evolution of the scatter and the mean 
Eddington ratio. 

3.2. Quasar Clustering 

With the model parameters constrained by the quasar LF, 
we are now able to make predictions for the clustering of 
quasars as a function of luminosity and redshift. Recall that 
our model is characterized by two parameters, the quasar life- 
time, tq, and the normalization of the Mbh-M/, relation, a. 
In the model, we assume that quasars are a random sample of 
the BHs in halos, and therefore tg has no effect on the cluster- 
ing of quasars. The clustering is therefore only sensitive to a, 
and this parameter is well-constrained at z < 4 (see Figure|4|. 
Moreover, a has an increasingly minor effect on the predicted 
clustering at higher redshifts. 
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Figure 6. The projected correlation function, Wp{R), vs. projected distance, 
R, at 5 redshifts chosen to be representative of the data. We include results 
fromUpss et al. ( 200l R09),^hite et al. 1 2012^ W12), and Shen et al. 1 2003 
S09), all of which are based on data from the Sloan Digital Sky Survey. At 
the highest redshift there is some tension between the model and data, but 
the eiTor bai's ai'e large and the simulation box is too small to provide model 
predictions at the largest scales. Future measurements of the clustering of 
both low and high redshift quasars will provide powerful constraints on the 
model. 

Figure |6] shows a comparison of our model and the data on 
the projected autocorrelation function, Wp{R), as a function 
of projected (comoving) distance, R, for a variety of redshifts 
chosen to illustrate the current constraints. We have computed 
the model correlation function by populating the halos drawn 
from an N-body simulatiorQ with BHs using the best-fitting 
relations derived above, and then calculating the clustering of 
BHs within the luminosity range of each observational sam- 
ple. This allows us to take into account the scale-dependent 

^ The simulation employed 2048^^ particles in a cubic box of side length 
1 Gpc with a force softening of 14kpc (comoving) and was run with the 
TreePM code of White 1 2002). Halos were found with a friends-of-friends 
algorithm (Davis et al. 1985) with a linking length of 0.168 times the mean 
inter-particle spacing. Spherical over-density masses were computed for 
each halo (including a correction for finite resolution). For the range of 
halo masses and redshifts of interest, masses defined via 180x the back- 
ground density are almos t identical to the 'virial' definition employed by 
iBehroozi, Wechsler & C onro^ f20O) 
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Figure 7. The large-scale bias predicted by our model as a function of lu- 
minosity for a number of redshifts. The relation is shallow at low luminosity 
due to the steepness of the Mbh relation at low mass (see Figure (8). 
The steepness of the relation at high luminosity depends on the scatter in the 
model, being less steep for more scatter. We have marked on the curves where 
the quasar number density is 5 X 10^'Mpc"', which con'esponds to of order 
100 quasar pairs within 20Mpc in a survey volume of lO'^Mpc^. To accu- 
rately measure the bias of objects at lower space densities one would need to 
resort to cross-correlations. 



bias and non-linearities, which are important on Mpc scales. 

We have populated the halos in the simulation at random, 
neglecting any properties of the halos apart from their mass 
(e.g., whether they have had a recent major merger). The 
majority of models assume that quasar activity occurs due to 
the major merger of two gas-rich galaxies, since this scenario 
provides the rapid and violent event needed to funnel fuel to 
the center of the galaxy (e .g. via the bars-within-bars instabil- 
ity; |Shlosman]etal] [1983) and feed the central engine while 
at the same time giving a connection between black hole fu- 
eling and th e growth of a spheroidal stellar component (e.g., 
iHopkins et al. 2008) . If black hole growth is feedback lim- 
ited, it is only a rapidly growing potential well that can host a 
rapidly accreting BH. 

Fortunately for such theories, the probability that a halo 
will undergo a major merger in a short redshift inter- 
val is only weakly dependent on the mass of the halo 
(Lacev&Cole 1993; Percival et al. 2003; Cohn & Whit 
2005;' Wetzel Cohn & White .2009,; .Fakhouri & Ma .200' ' 
Hopki ns et al. 1 1201 OH) , i.e., the mass function of such ha- 
los is almost proportional to the mass function of the par- 
ent population. Similarly the clustering properties of recently 
merged halos are similar to a random sample of the popu- 
lation with the same mass d istribution ( Percival et al.. ,20031 ; 
I Wetzel. Cohn & White! l2009h . Thus in any redshift interval 
the fraction of halos of mass M/, that undergo a quasar event 
is almost independent of M/, and z and can be regarded as a 
random selection. 

The agreement between the data and the model is excel- 
lent at z < 3, especially considering that the model was only 
tuned to the quasar LF. The model under-predicts the cluster- 
ing at z ^ 3.7, although the errors on the data are large. This 
model prediction is quite robust: the Mbh-^sI relation at 
high redshift becomes very steep (see Figure [8] below), and 
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Figure 8. The typical black hole mass in the central galaxy of a halo of mass 
Ml,, vs. Ml,, for a number of redshifts (corresponding to the redshifts shown in 
Figure[3). The typical BH mass con'esponding to a fixed M;, increases with 
z, as expected. Note the significant curvature in the relation, which arises 
due to our assumption that galaxy properties regulate the size of black holes 
and the well-known inefficiencies of galaxy formation in high and low halo 
masses. Since the .v-axis controls the clustering amplitude while the v-axis 
controls the observed luminosity, the curvature of this relation has im porta nt 
implications for the luminosity dependence of quasar clustering (see q3.2) . 

SO even a significant change in a or 77 changes the clustering 
only modestly. Similarly, changes in the assumed Lq-Mbh 
scatter within the range 0.3-0.6 dex does not significantly al- 
ter the predicted clustering. This occurs because a change in 
scatter induces a change in a that happens to leave the clus- 
tering essentially unchanged. Future constraints on the clus- 
tering of high-redshift quasars will place strong constraints on 
this model, as discussed further in i i4.1l and may indicate that 
some of our model assumptions break down as we approach 
an era of rapid BH growth at high z- 

Observationally, it has proven very difficult to measure a 
dependence of clustering strength on quasar luminosity (see 
e.g., IShen e t al. 200^ for a recent example), in part be- 
cause the significant scatter between quasar luminosity and 
halo mass will dilute any intrinsic relation between cluster- 
ing strength and luminosity. We address this issue in Figure 
|7] where we plot the large-scale bias as a function of lumi- 
nosity and redshift. Here the model bias was computed via 
the relation betwe en bias, halo mass, and cosmology from 
[Tinker et alldMoh . 

We find a very shallow relation between bias and quasar 
luminosity below M,(z = 2) -26. In our model this occurs 
for three reasons: (1) the intrinsic relation between bias and 
halo mass is very shallow below the characteristic halo mass, 
which at z is IQ^^Mq; (2) the Mbh-M,, relation be- 
comes very steep at low mass, implying that a large range in 
quasar luminosities maps into a small range in halo masses; 
(3) scatter in the Mgai-M/,, M^n-Ma-^x, and Lq-Mqu relations 
acts to dilute the strong clustering in high mass halos. The de- 
gree of luminosity dependence (as well as the absolute value 
of the bias) is sensitive to the scatter in the Lq-M/, relation, 
with more scatter leading to less L-dependence. This weak 
luminosity-dependent clustering is also predicted in the mod- 
els of .Croton (.2009.) . . Shea C2009I) . and lHopkins et al.. C2008I) . 
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Figure 9. Relation between halo mass and redshift for quasars of a fixed lu- 
minosity. At low redshift the range of halo masses hosting quasars is very 
broad, but the distribution naiTows substantially at high redshift. This is sim- 
ply a recasting of the relations shown in Figure[8] 

Figure|7]demonstrates that we expect significant luminosity 
dependent quasar bias only for very luminous quasars. How- 
ever, measuring the autocorrelation function of such luminous 
quasars is made difficult by their low space densities, which 
can be illustrated as follows. The error on the bias in the 
high-L regime is dominated by counting statistics. The num- 
ber of pairs within e.g., 20Mpc is (1 /2)fiQ [1+^20] VsurveyV2o 
where V20 = (47r/3)(20Mpc)-', Vsurvey is the survey volume, 
HQ is the quasar space density, and ^ is the volume aver- 
age correlation function. For ^(r) = (ro / r)^ we have ^ = 3^, 
and ro ^ 10-20/!~'Mpc so we expect ^ ^ 0{1). One hun- 
dred pairs within 20 Mpc would return an error on the bias of 
^ 10%, and for a fiducial survey volume of 10"^Mpc^^, this 
corresponds to a quasar number density of ~ 5 x 10~^Mpc~^. 
The luminosity corresponding to this number density at each 
redshift is marked by a solid symbol along the b{L) relation 
in Figure |7] In order to probe the bias for quasars at higher 
luminosities it will be necessary to resort to cross-correlation 
techniques, which allow estimates of the bias of objects with 
extremely low space density. An appealing method would be 
to cross-correlate existing spectroscopic samples of quasars 
with samples of galaxies or lower luminosity quasars selected 
from deeper photometry in upcoming surveys such as DES, 
Pan-STARRS, SUMIRE and LSST. 

Finally, we emphasize that there are several model param- 
eters that were assumed herein that affect the clustering of 
quasars. In particular, the scatter between Lq, Mgai and M/, 
has a significant effect on the luminosity-dependence of bias. 
As mentioned previously, the Eddington ratio, 77, is perfectly 
degenerate with a, and so with the freedom to vary the nor- 
malization in the MaH-A^gai relation, the clustering will be 
insensitive to rj. Seen another way, constraints on a from 
observations will in turn allow a constraint on 7/ within the 
context of our model. 
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4. DISCUSSION 
4.1. Implications 

The success of our model in reproducing the basic demo- 
graphics of quasars allows us to consider several implications 
that follow naturally within our framework. 

In Figure [8] we show the best-fit model Mi,-Mbh relations 
from z = 0.5 to z = 3.75 (the model at z = 4.75 is highly under- 
constrained and so is not plotted). As discussed above, the 
quasar LF places very weak constraints on the model rela- 
tions at log(M/,/MQ) > 13.5, and so one should interpret the 
model relations in Figure[8]with this in mind. It is also worth 
pointing out that while the model formally allows for the ex- 
istence of extremely massive BHs with Mbh > 10'"Mq re- 
siding within moderately massive halos, at high redshift such 
halos are very rare. For example, at z = 4.75 one expects only 
of order one halo with log(M/,/Mo) >13 per 10^ Mpc"'. 

In the model, lower mass black holes are growing to lower 
redshift faster than higher mass black holes (this is some- 
times referred to as downsizing). In fact, a strict inter- 
pretation of the results has the disconcerting property that 
the BHs in the most massive galaxies lose mass below z ~ 
2. This arises directly from two observational facts: (1) 
MBH/Maai declines with time as ( 1 +zf (M cLure et alJl2006l: 
Ir^rgett. Dunlop & McLurell2012h and (2) the most massive 
galaxies grow very little between z ~ 2 and the present, at 
least within the inner several kpc where high-redshift galaxies 
are typica lly photomet e red (e.g ., Perez-Gonzalez et al. 2008; 
iMarchesini et alJl2009l: Ivan Dokkum et all 12010). There are 
several caveats worth noting in this context. First, the slope 
of the Mbh - Afgai relation is not well-known at high Mgai, ow- 
ing largely to the rarity of such objects. Second, the diffi- 
culty in measuring the total light from galaxies implies that 
non-trivial systematic uncertainties still plague conclusions 
regarding the evolution of the most massive galaxies from 
< z < 2. Finally, as demonstrated in Figure [5] the model 
is essentially unconstrained at the highest masses, and so any 
conclusions regarding evolution at the high-mass end must be 
treated with caution. Evolution in Eddington ratios (advo- 
cated by Shen & Kelly 2012), curvature in the scaling rela- 
tions or strong evolution in bulge-to-total mass ratio could all 
play a role in alleviating the tension. Future observations of 
the evolution of the MeH-Mgai relation at high mass will be 
necessary to make further progress. 

Figure |9] shows the evolution of the halo mass for quasars 
of fixed luminosity. The trend of lower M/, at higher z was al- 
ready apparent in Figure[8] Figure|9]also emphasizes how the 
range of halo masses for a fixed luminosity range narrows to- 
wards higher z. This effect is in the opposite sense to models 
whic h tie the lum inosity of quasars directly to halo properties 
(e.g. ICrotonI 12009 ). Our model is able to reproduce the ob- 
served L-independent clustering at low z because the run of 
bias with halo mass also becomes shallower at low z for the 
halo masses of interest. 

Figure [TO] shows the evolution of the quasar luminosity 
function as predicted by our model once the parameters are fit 
to the available data. The evolution of the LF is driven by evo- 
lution in the Mbh - Afgai and Mgai -M/, relations and the evolu- 
tion of the halo mass function (the Lq -Mbh relation does not 
evolve in our fiducial model). The break in the model quasar 
LF arises primarily due to the shape of the Mgai -M/, relation, 
and thus quasars live in halos near the peak of that relation. 
Ml, ^ IO'^Mq. The peak of the Mg^\-Mi, relation changes 
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Figure 10. The best-fit quasar luminosity function for a variety of redshifts. 
Notice the rapid drop in the LF at high luminosity at z < 2, which, in our 
model, is largely responsible for the rapid drop in the cosmic quasar number 
density at late times. The break in the quasar LF scales approximately as 
(l+z)", see the text for details. 

very little with redshift (Behr oozi. Wechsler & Con rov''2012l 
e.g.,), so that at fixed M/, there is little change of Mgai with 
z. However the luminosity of the break can evolve due to 
a combination of evolution in the Mbh -Mgai relation or the 
Eddington ratio. In our fiducial model rj is constant and 
Mbh fx (1 +z)^ at fixed Mgai and so the break in the luminosity 
function scales as (1 +z)^. The faint-end slope of the model 
LF does not vary significantly, in good agreement with the 
data, and the overall normalization changes only modestly. 
The major departure from pure luminosity evolution is the 
change in the slope of the bright end. The bright-end slope 
appears shallower at higher z both because the data are prob- 
ing closer to the (brighter) break of the LF and because the 
Mbh - Ml, relation becomes steeper at higher mass and red- 
shift. We also note that the bright end of the model LF is 
strongly suppressed at z < 1 .5, and it is this suppression that is 
responsible for much of the drop in the quasar number density 
to lower redshift. The drop is a consequence of the shallowing 
of the Mbh- Ml, relation at high mass, which is in turn driven 
by the very slow growth of massive galaxies at low redshift. 

In fact, the model naturally reproduces the global rise and 
fall of the quasar number density over the interval 0.5 < z < 
4.75. This follows simply from the evolution in the Mgai-M/, 
and Mbh -Mgai relations and the halo mass function; it does 
not require strong evolution in Iq at low z. Specifically we 
do not invoke a decline in the cold gas fraction nor a decline 
in the major merger rate at z < 2 in order to reproduce the 
observed decline in the abundance of quasars. While these 
physical processes may ultimately be responsible for shaping 
the evolving relations between Mbh, A^gai and M;,, they do not 
appear explicitly in the model. 

Our model suggests a very different picture of how quasars 
inhabit massive halos compared to previous work. Rather 
than having a preferred halo mass scale (around IO'^Mq) for 
quasar activity we have actively accreting black holes in a 
broad range of galaxy and halo masses. The apparent pref- 
erence for quasars to live in halos of lO'^M© arises from 
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the shape of the Mgai-M;, relation, which reflects the wefl 
known fact that galaxy formation is most efficient in halos 
near IO'^Mq, along with the shape of the halo mass function. 

Due to its simplicity the model predicts the clustering 
of any population of quasars once the two model parame- 
ters are fixed (e.g., by the observed LF). Variation in the 
Lq - Ml, scatter or Mbh - A^eai slope do not strongly af- 
fect the predicted clustering, meaning that our model makes 
an essentially parameter-free prediction of the clustering of 
quasars as a function of luminosity and redshift. Over- 
all the agreement between the predicted clustering and the 
observations is good, though there is a tendency for the 
model to slightly underpredict the observations and there 
is some tension at the highest redshifts. This tension has 
been noted before - the very high amplitude of clustering 
measured at z 4, in combination with the abundance, re- 
quires quasars to have a duty cycle approaching uni ty and 
almost no scatter in Lq at fixed Mh (, White. Martini & CohnI 
120081 iShankar. Weinberg & Shenll2010l) . If these measure- 
ments can be strengthened, possibly by cross-correlation of 
existing spectroscopic quasar samples with deeper photomet- 
ric quasar or galaxy samples, then it will indicate that one of 
our assumptions is breaking down as we approach the era of 
rapid black hole growth in the early Universe. 

We make no assumption about what triggers quasar activ- 
ity, whether it be a major merger of two gas rich galaxies, a 
secular instability in a disk, or a critical halo mass. In gen- 
eral it is quite difficult to translate abundance and clustering 
measurements into constraints on the underlying mechanisms 
that trigger quasar activity. We can gain some insight by the 
fact that our duty cycle, or quasar lifetime, is relatively in- 
dependent of redshift with a tendency to fall towards higher 
redshifts rather than rise. If quasars are visible for a fixed, 
but short, time and are triggered by mergers then we expect 
fg to scale with the merger rate (c.f. Carlberg 1990). The 
mer ger rate for halos, per halo, per unit redshift is relatively 
flat dLacev& Cole"! 993^, 'Percival et al."2003^, 'Cohn & Whitd 
20051: IWetz el, Cohn & White 2009; Fakhouri & Ma 200i 



Hopkins et al.l i2010tal) . so if we can naively translate halo 



mergers into galaxy mergers we expect a rate (per unit time) 
scaling as ( 1 + z)H(z) oc ( 1 + z)^/^ for z ^ 1 . If quasars are vis- 
ible for a constant interval after each merger then fg cx 1+z, 
which is not in good agreement with our best-fit relation. Of 
course, galaxy merger rates can differ from halo merger rates. 
A recent analysis by Hopkins et al. (2010a) suggests a rate per 
unit time scaling as (1 +z)'-^"^ '^, which would lead to slower 
evolution in tg, as we observe. Such agreement is not conclu- 
sive however, and we cannot rule out secular processes or a 
time-varying combination of multiple triggers. 

4.2. Comparison to Previous Work 

The success of our model in explaining the basic demo- 
graphics of quasars with relatively few, smoothly varying 
inputs goes a long way to explaining the manner in which 
forward modehng of the quasar population can succeed with 
relatively Uttle fine tuning Both s er ni-analytic models (e.g., 
Cattaneo. Haehnelt & Rees* '1999^, 'Kauffmann & Haehnelt' 
200a I200I Volonteri, Haardt & Madau 2003; 
Bromlev, Somerville & Fabian 2004; Granato et al. 2004; 
Croton et al., ,2006; Monaco, Fontanot & Taff oni> , 2007 1 



i 



Malbon etal.]|2b 07; Bonoli et al. 2009; Fanidakis et al.ll201 
Hirschmann et al. 2012) and hydrodynamic simulations (e.g., 
Siiacki et al...2007; ,De Graf et al...201 1.) adjust their subgrid 



models to ensure a reasonable match to the M„^\-Mh relation 
over a broad redshift range, thus ensuring that galaxies 
populate halos in approximately the correct manner. All of 
the models introduce a Mbh - gai relation through either 
or a combination of common feeding mechanisms and 
feedback-limited BH growth. As we have shown, with these 
two ingredients even simple lightcurve models are sufficient 
to match the basic demographics of quasars over a broad 
range of luminosity and redshift. A good match to the data 
can be found for a wide range of scatter in M^n- M„s.h or 
evolution in the scatter. Conversely, if a model has difficulties 
reproducing the stellar mass function and its evolution then it 
will need to incorporate mass-dependent quasar physics that 
counteracts this deficiency in order to match the observed 
quasar properties. 

By contrast, models that tie black hole properties directly 
to the underlying halo population need to introduce more 
complexity in order to reproduce the observed prop e rties of 
quasars. Re cent e xamples include iLidz et alj (120061) . ICrotonI 
(I2b09i) . and lShenI (12009). who all need to include mass- and 
redshift-dependent duty cycles to explain the shape and evo- 
lution of the quasar luminosity function. While our model and 
theirs can produce qualitatively similar fits to the basic data, 
the explanations for the observed behaviors differ. One of the 
most basic differences is the range of halos that host active 
quasars, and its evolution (discussed above). This in turn af- 
fects how each model explains the evolution of the quasar LF 
and the luminosity-independence of quasar clustering. 

Conventional wisdom is that the quasar duty cycle is re- 
quired by the data to be a (strong) fu nction of luminosity (e.g. 
Adelberger & Steidel 2005; Hopki ns et all I200I iLidzetal] 
l2006i: ■Croton„2009; .Shen.,2009 ). In our model this is not the 
case. There are two major reasons for this. The first is that 
we obtain a flattening of the b(L) relation from the steepness 
of the Lq - Mh relation at low Lq and second is the intrin- 
sic scatteJ3 in that relation. Thus our model is n ot a "light 
bulb" model in the sense o f iHopkins et al] (l2005h ; lOdz et alj 
(I2006h . who reserve that term for a model in which there is 
no scatter in Lq-M/,. However scatter in the Lq-M/, rela- 
tion is expected, due to the observed scatter in Mbh - Afgai 
and variation in Eddington ratios if from no other source; for 
this reason we refer to our model as a "scattered" light bulb 
model. This expected level of scatter is enough to make b(L) 
flat until extremely high L or correspondingly low ng (a sim- 
ilar behavior is seen in the model of Croton 2009, which is 
also not strictly a light bulb model in the above sense). For 
this reason we are able to obtain a model in which both the 
quasar lifetime and the quasar clustering are independent of 
L 

lAird et aP (120 12h studied X-ray selected active galactic nu- 
clei (AGN) as a function of galaxy mass at z ^ 0.6 and found 
no preference for AGN to be found in galaxies of a particu- 
lar mass at fixed Eddington ratio, even for ratios as high as 
77 > 0.1. Their results suggest a duty cycle that does not de- 
pend strongly on galaxy mass, in excellent agreement with 
our results. 

Finally, the apparent preference for quasars to live in ha- 
los of 1O'^M0, which has been noted by many authors, arises 
in our model from the shape of the Mgai -M/, relation which 
reflects the well-known fact that galaxy formation is most ef- 

^ This scatter may arise due to time-dependent processes, i.e. a high Lq ob- 
ject at the time of observation is not required to have always been or continue 
to be high Lq. 
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ficient in halos of IO'^Mq. Within the context of our model 
this cannot be taken as evidence for a merger driven origin 
to quasar activity, despite the fact that it is close to the small 
group scale where mergers may be more efficient, because it 
is not believed that the peak of the Mg^i-M/, relation at this 
scale is related to mergers. 

4.3. Mock Catalogs 

While our intent has been to understand the quasar phe- 
nomenon, the model can also be used for the creation of mock 
catalogs from N-body simulations. The simplicity of the 
model makes it easy to rapidly generate redshift-dependent 
quasar populations that have the correct luminosity function 
and clustering, given halo catalogs at the redshifts of interest. 
The steps for creating such a catalog are straightforward: 

1. Adopt the redshift-dependent Mgai - M/, relation from 
iBehroozi. Wechsler & Conrovl(l2012l) . including scatter 
inMgai at fixed M/,. 

2. Use the Mgai-Men relation from Equation[T]to assign 
BHs to galaxies, including 0.3 dex of scatter in Mbh at 
fixed Mgai. As our be st-fit normalization is close to the 
iMcLure et al.l (120061) relation, especially at z < 3, we 
suggest fixing a = -3.1 in this equation. 

3. Randomly turn a fraction, of the BHs into active 
quasars. As evident from Figure HI the quasar lifetime 
is approximately constant at 3 x 10^ yr at z < 3; we 
therefore advocate fixing Iq to this value. One then de- 
termines the duty cycle via foniz) = tg/tHiz). 

4. For the active BHs, convert Mbh intoLg using Equation 
121 with an additional 0.3 dex of scatter in Lg at fixed 

When simulations are populated with quasars in this way, 
the mock quasar LF and clustering will agree with all existing 
LF and clustering data at z < 3. In order to produce mock 
catalogs at higher redshifts one will need to include a drop in 
tg as shown in Figure ID Such mock catalogs should prove 
useful in the context of ongoing and future planned surveys 
such as BOSS, bigBOSS, DES, Pan-STARRS, SUMIRE and 
LSST. 

5. SUMMARY 

We have presented a simple model for quasars with the aim 
of understanding to what extent their demographics arise nat- 
urally from what is known about the evolution of galaxies 
along with plausible assumptions about how black holes in- 
habit them. The key feature of the model is that the properties 
of black holes are set by those of their host galaxies rather than 
their host halos (see also iWhite et al.l 120121) . In the model, 
BH mass is linearly related to galaxy mass and BHs shine at 
a fixed fraction of the Eddington luminosity during accretion 
episodes. Galaxies are related to dark matter halos via em- 
pirically constrained relations ( Beh roozi. Wechsler & Conrovl 
I2OI2I) . The model has only two free parameters at each red- 
shift, the normalization of the MfiH-Mgai relation and the duty 
cycle, both of which are tightly constrained by observations of 
the quasar LF. The model provides an excellent fit to the LF 
data for 0.5 < z < 5 and reproduces the observed clustering 
at intermediate redshifts with no additional adjustable param- 
eters. 



The best-fit model parameters have Msu/Mg^] cx (1 + z)^, as 
observed, and a quasar lifetime of approximately 10^ yr in- 
dependent of epoch. This may be expected if the growth of 
the galaxy during a quasar event only allows ^ 1 e-folding of 
black hole growth before feedback halts quasar activity. 

There are several implications of our model, which we now 
summarize: 

• Actively accreting BHs are equally likely to exist in 
galaxies, and dark matter halos, over a wide range in 
masses. The BHs in halos more massive than IO'^'^Mq 
contribute very little to the observed quasar LF at any 
redshift due to their rarity. The quasar LF therefore 
places weak constraints on the quasar duty cycle in 
massive halos. 

• The break in the quasar LF is a reflection of the peak in 
the Mgai-M/, relation at M/, ^ lO'^M© and the observed 
evolution of the LF primarily reflects the (1 +z)^ scaling 
of Lg /Mgai and the change in shape of the Mgai - M/, 
relation. The bright-end slope of the quasar LF appears 
shallower at high z both because the data are probing 
closer to the (brighter) break in the LF and because the 
Mbh-Mii relation becomes steeper at higher mass and 
redshift. 

• Our model naturally reproduces the global rise and fall 
of the quasar number density over the interval 0.5 < z < 
5. This follows simply from the evolution in the Lg- 
Mk relation and does not require strong evolution in the 
quasar lifetime. The bright end of the model quasar LF 
is strongly suppressed at z < 1.5, due to the slow growth 
of massive galaxies, and this is responsible for much of 
the drop in quasar number density to low redshift. 

• The apparent preference for quasars to live in halos of 
1O'^M0 arises from the shape of the Mgai-M/, relation, 
which reflects the well-known fact that galaxy forma- 
tion is most efficient near IO'^Mq, in conjunction with 
the steepness of the halo mass function at high mass. 

• There is some tension between our model and the am- 
plitude of clustering observed at z ~ 4; the latter, taken 
at face value, suggests that quasars have a duty cycle 
approaching unity and almost no scatter in Lg-M/,. Fu- 
ture measurements in this redshift range will be crucial 
tests of the model. 

• The nearly constant inferred quasar lifetimes as a func- 
tion of luminosity and redshift should provide valuable 
constraints on the triggering mechanisms for quasars. 

Measurements of quasar demographics at higher redshifts 
and lower luminosities will help to further constrain and test 
our model. In particular, stronger constraints on the quasar 
LF at z > 4, quasar clustering as a function of luminosity and 
redshift, and the Mbh - Mgai relation as a function of redshift, 
will provide very strong constraints on the model parameters. 
Moreover, with such observational constraints in hand, we 
will be able to directly constrain the mean Eddington ratio 
as a function of redshift and the scatter as a function of red- 
shift, providing further insight into the link between quasars, 
galaxies, and dark matter halos. 
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